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With the standard deviation for the logarithm of the re-scaled range {\F{t+T) — F{t)\) of simulated 
fractal Brownian motions F{t) given in a previous paper [I[, the method of least squares is adopted to 
determine the slope, S, and intercept, I, of the log({|E(t-|-T) —-F(t)|)) vs log(r) plot to investigate the 
limitation of this procedure. It is found that the reduced of the htting decreases with the increase 
of the Hurst index, FI (the expectation value of S), which may be attributed to the correlation among 
the re-scaled ranges. Similarly, it is found that the errors of the fitting parameters S and / are usually 
smaller than their corresponding standard deviations. These results show the limitation of using 
the simple least square method to determine the dimension of a fractal time series. Nevertheless, 
they may be used to reinterpret the htting results of the least square method to determine the 
dimension of fractal Brownian motions more self-consistently. The currency exchange rate between 
Euro and Dollar is used as an example to demonstrate this procedure and a fractal dimension of 
1.511 is obtained for spans greater than 30 transactions. 
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I. INTRODUCTION 

The properties of fractional Brownian motions (fBms) have been investigated by researchers in different fields, e.g. 
statistics, hydrology, biology, finance, public transportation and so on, which has helped us to better understand many 
complex time series observed in nature The Hurst exponent H {0 < H < 1) is the most important parameter 

characterizing any given time series F{t), where t represents the time steps, and the fractal dimension D is determined 
via the relation D = 2 — H. The Hurst exponent H is defined with the following expression [I|: 

log((|F(t + r)-F(t)|))ociLlog(r), (1) 

where '()' represents averaging over t. Using the Lowen method for 0 < < 0.5 and the circulant embedding 

method for 0.5 < H < 1, Qiao & Liu (2013) carried out extensive simulations of fBms with different H to obtain 
the standard deviation of the re-scaled range (| A Fi-|) = (|F(t -|- r) — F{t)\) for different sampling methods. This 
paper extends this study to investigate the limitation of applying the commonly used least square (LS) method to 
the log((|E'(t -b t) — F{t)\)) vs log(r) plot for the evaluation of H. 

This paper is organized as follows. With the sampling method 4 of Qiao & Liu (2013), the standard LS method 
with weight is applied to the log((|F(t -I- r) — F(t)|)) vs log(r) plot for a few fBms with different FI in Section 2, and 
the results of the fitting are analyzed. In Section 3, a similar work is carried out for the third sampling method of 
Qiao & Liu (2013). Einally, we give our conclusion and a brief discussion in Section 4. 


II. THE LS FITTING FOR ONE SAMPLING METHOD 


In the following, we use discrete-time fBms F{t) to obtain the dependence of the rescaled range (|AFAt|) = 
(|F(t + At) — F{t)\) on the time span At. The sampling method corresponds to the case 4 of Qiao, B.Q and Liu, 
S.M, (2013), which has the lowest standard deviation for the rescaled range and can be expressed as 


(IAFa*!) 


Int(Ar^/Ai)-l 

^ [|E'(At-b At * f) — F(At =1= f)|] 
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where N = 250 and At is a positive integer. According to our previous work [ij, the standard deviation of log((| AEAt|)) 
is given by: 


cr = A-Hog(e)(7r/2 - l)A2(At)A2 for 0 < iL < 0.5 
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and 


o- = 7V-%g(e)(7r/2-l)i/2(^^)/3/2 0.5 < i? < 1.0 

where [3 = 2 — 2H + 0.26^^^ — [i2H — 1.5)^ + 0.01]^/^. It should be noted that all of the measurement errors are only 
in the log((|AFAt|)) variable and the values of At are all exact without errors. Then the standard LS method with 
weight (e.g. chapter 15 of Press et al. 0) can be applied to the log((|F(t + At) — F{t)\)) vs log(At) plot to 
obtain the slope S, the intercept /, and their corresponding errors Se, A, respectively. Figure [T] shows these results 
for several values of H. 



FIG. 1: Left: examples of a few fBms with different values of H. Right: dependence of the re-scaled range (|AFAt|) on the 
time span At for H = 0.8 (magenta), H = 0.5 (green), H = 0.2 (red). The three black solid lines give the best fit obtained 
with the standard LS method. The fitting parameters and errors are indicated in the figure. The normalized residuals r of 
log((|AFAt|)) for different H values are shown in the three panels below. 



FIG. 2: Distribution of the normalized residual r for H = 0.8, 0.5, 0.2 shown in FIG. 1 and the Gaussian fitting is indicated by 
the red solid line with the fitting parameters also given in the figure. 

We first note that the normalized residuals r = ((|AFAt|) — {|AFAt|))/cr, where (|AFAt|) corresponds to the best 
fit model, do not follow the standard Gaussian distribution as expected. Figure [2] shows the distribution of these 
residuals and the corresponding fit with a Gaussian model. The width of the Gaussian fit is less than I for the three 
values of H studied here. The expectation values of S and I are H and log(2/7r)/2(~ —0.0981), respectively 0 . 
Although the best-fit values are close to these expectations, considering the small errors of these fitting parameters, 
the deviation of the best-fit values from their expectation appears to be significant. 

To better understand these issues, for 19 different uniformly sampled values of H between 0 and 1, the simulation 
and fitting above are repeated 2000 times. Figure [3] shows the dependence of the square-root of the reduced on H 
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FIG. 3: The H dependence of the square-root of the reduced the Gaussian width of the normalized residual, Or, for 20 

runs randomly selected from 2000 runs. The solid squares indicate the mean value of the square-root of the reduced x? for the 
2000 rims. 



FIG. 4: Left: the distribution of the best-fit slope for 2000 runs with H = 0.8 and the Gaussian fit. Right: same as the left 
panel but for the intercept. 


for 20 randomly selected simulations from the 2000 runs for each value of H. The stars show the corresponding widths 
of the Gaussian fits to the normalized residual, Ur, as shown in FIG. 2, which are consistent with the value of the 
reduced ■ We note that both the reduced and the width of the distribution of the normalized residual decreases 
with the increase of H. Figure U show the distributions of the best-fit intercept and slope and their corresponding 
Gaussian fit for H = 0.8. 

The distributions of S and I do not necessarily follow Gaussian distributions, especially for values of H close to 0 
or 1. In the following, we use the standard deviation of S (denoted with as) and I (aj) to characterize the widths 
of these distributions. Figure [S] shows the dependence of /e/c/ and S^/as on H. The widths as and ct/ are usually 
greater their errors Se and le shown in FIG. I, respectively. We note that Se and /g are determined by the standard 
deviation cr and therefore do not change for the 2000 simulations and fittings with a given value of H^. Figure [S] also 
shows the correlation coefficient R between Ig and Sg-, which corresponds to the non-diagonal term of the covariance 
matrix of the fitting parameters S and I. The fact that R is close to —1 implies that the two fitting parameters are 
not independent. The anomalous behavior of Igjaj and Sgjas for H > 0.7 may be related to the increase of the 
correlation between Ig and Sg as shown with R. 

The above results clearly show inconsistency in the evaluation of H with the LS method applying to the log((|F(t-|- 
At) — F{t)\)) vs log(At) plot. The strong correlation between the errors of the fitting parameters suggests that 
correlations among log((|F(t -|- At) — F(t)|)) may play a role. For the simple LS method to be applicable, each 
rescaled range (lAFAtl) needs to be independent of the others, which does not appear to be the case in our analysis 
since the rescaled ranges are derived from a given fBms. The upper panels of Figure [5] show the correlation between 
log((|AF 2 |)) and log((|AF 3 |)) for several values of H. The lower panel shows the dependence of the correlation 
coefficient 7 on iJ. The correlation increases with the increase of H, reminiscent of the dependence of the reduced 
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FIG. 5: Dependence of the 7e/o-/, Se/(^s and R on H for 2000 runs of each H. 







H 


FIG. 6 : Upper: correlation between log((|F’(t+2)—F(t)|)) and log({|F(t+3)—F’(t)|)) for 2000 runs of fBms with H = 0.8, 0.5, 0.2. 
The corresponding correlation coefficients are indicated on the figure. Lower: dependence of the correlation coefficient 7 in the 
upper panel on H. 


of the fitting on H shown in Figure [31 


III. RESULTS FOR ANOTHER SAMPLING METHOD 

Qiao and Liu (2013) showed that the standard deviation of the rescale range has relatively simple expression for 
the following sampling method (case 3 of i): 
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The standard deviation of log((|AFAt|)) is given by 

cr = iV“i/2log(e)(7r/2 - for 0<H < 0.5, 

cr = N-^^^\ogie)iTT/2 - f)i/2 for 0.5 < H < 1.0, 

which are larger than those for the sampling method in Section 2. Figures 7-9 show the corresponding results, which 
are very similar to those in Section 2. The most notable difference appears in R and 7. With the current sampling 
method, the (anti-)correlation between the errors of the fitting parameters Se and R is even stronger, while the 
correlation between log((|AF 2 |)) and log((|AF 3 |)) decreases slightly. 


Fractal Brownian Motion 
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FIG. 7: Same as FIG. [T] but for the sampling method given by equation |3] 




FIG. 8: Left: same as Fig[3]but for the sampling method given by equation[3l Right: same as FIG.[5]for the sampling method 
given by equation [3] 


IV. CURRENCY EXCHANGE RATE BETWEEN EURO AND DOLLAR 

The results above show the limitation of the least square method in the evaluation of the dimension of fractal time 
series. However, we notice that the normalized residual of the LS fit shown in Figures. [1] and 0 does not appear to 
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FIG. 9: Same as FIG. [6] but for the sampling method given by equation |3] 


depend on the time span, which also suggests the soundness of this procedure. One therefore may reinterpret the 
fitting parameters to obtain more self-consistent results. In this section, we use the currency exchange rate between 
euro and dollar in May and June of 2008 as an example to demonstrate how to obtain the dimension of fractal time 
series self-consistently. 

The top-left panel of Figure fTOl shows the exchange rate as a function of the number of transaction t. It covers a 
duration from May 1 to June 30. The market was closed during the weekends. In total we have 43 days with extensive 
transactions. Since the sampling method corresponds to the case 4 of Qiao, B.Q and Liu, S.M, (2013) has the lowest 
standard deviation for the rescaled range, we will only consider this case in the following analysis. 



FIG. 10: Left: Exchange rate between Euro and Dollar between May 1 and June 30 and its Hurst index analysis. The red 
open circles show the dependence of the square root of the reduced (indicated by the scale on the right) on the lower bound 
of the fitting range. The pink line indicates the best fit result with a lower bound of 30; Right: The Hurst index of exchange 
rate between Euro and Dollar for 43 days with extensive transactions. 

The bottom-left panel of Figure [TO] shows the dependence of the rescaled exchange rate range on the span of the 
transactions At. Because the error bars are very small, they are almost invisible in the figure. It is evident that 
the rescaled exchange rate range has some curvature below a few tens of transactions. Within ~ 4 transactions, the 
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rescaled range does not appear to vary with the transaction span, implying complete randomness similar to white 
noises. However, above a transaction span of ^ 30, the rescale range can be fitted linearly with a slope of 0.489±0.002, 
corresponding to a fractal dimension of 1.511. The reduced = 0.48, which is consistent with results shown in Figure 

H 

To further demonstrate the consistency of the above procedure, we analyze the data for the 43 days separately. 
The top-right panel of Figure [TO] shows Hurst indexes, where the errors have been corrected using results in Figure 
[SJ The bottom-right panel of Figure [TU] shows the normalized residual, which is consistent with a standard Gaussian 
distribution. 


V. CONCLUSION AND DISCUSSION 

The scaling of the rescaled range (|AFV|) with the time span r plays an important role in quantitative analysis 
of fractal time series F{t). The Hurst exponent is defined as the exponent of a power-law scaling. However, the 
consistency of methods of the Hurst exponent measurement has not been explored. One of the commonly used 
method is to apply the lease square method to the log({|AF.r|)) vs log(r) plot. Using the mathematically well-defined 
fractal Brownian motions as a standard model, we investigate the limitation of the LS method. 

It is shown that, due to the correlation among the rescale range for the sampling methods adopted in the paper, 
the LS method does not give self-consistent results. The two sampling methods have the most comprehensive usage 
of the available data and relatively simple expressions for the standard deviation of the rescale range. Other sampling 
methods are relatively less practical and less efficient. In principle, one may remedy the above procedure by taking 
into account the correlation among the rescale range into account in the fitting process. In that case, instead of fitting 
a simple function log((|AFVl)) vs log(T), one also needs to know the covariance matrix of log((|AF.r|)) in advance. 
These calculations are above the scope of the current investigation. 

The fact that the normalized residual of the LS fit does not appear to depend on the time span also suggests 
the soundness of this procedure. The LS method may still be adopted in practical application by modifying the 
interpretation of the fitting result slightly with the results in this paper. First the reduced should have a best fit 
value given in figures [3] and [5] instead of I in the simplest LS method. Second the errors of the fitting parameters need 
to be corrected with the results shown in figures [5] and ID 
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